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ABSTRACT 


Section 208 of the Federal Water Pollution 
Control Act Amendments of 1972 provided the op- 
portunity and funding to fight -water pollution through 
the use of regional water quality planning* A 
common requirement of the 208 program is to de- 
velop a capability of predicting -water quality in the . 
rivers and lakes resulting from existing and po- 
tential land-use policies* To achieve this capability, 
the Ohio -Kentucky-Indiana (OKI) Regional Council 
of Governments is developing a deterministic model 
capable of predicting sediment and nutrient flow 
into the waterways. An essential input to OKI's 
model is an accurate map of land use within the 
watersheds. This information was obtained by OKI 
through the machine processing of LANDSAT-1 
digital tapes. Computer tabulations were generated 
to obtain area covered by each of 16 land-use cate- 
gories within 225 drainage areas* The 16 cate- 
gories were merged into 10 categories and mapped 
at a scale of 1 inch = 5, 000 feet with detail to 0.44 
hectares {1. 1 acres) for the 2,700-square mile 
region. The map products and data were produced 
within a period of less, than 9 0 days at a co st of 
$20, 000, a significant improvement in dollars and 
time over conventional mapping techniques . 


BACKGROUND 

State and federal agencies are becoming in- 
creasingly alarmed over the lo s s in .. water quality 
in many of our public lakes and rivers. Much of 
this loss is a direct result of pollution generated 
by man and the increased sediment and nutrient 
runoffs into the rivers and lakes resulting from; 
urbanization in the watersheds* It is now realized 
that our water resources are not inexhaustible and 
that land development in the watersheds must be 
planned if the conflict between utilization of our 
wate r re s our ces and maintenance of the quality, of 
our lives is to be resolved. 


To provide for this needed planning, the 
Federal Water Pollution Control Act Amendments 
of 1972 created several programs to fight water 
pollution* Under provisions of Section 208 of that 
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act (EPA, 1974), regional councils of government 
such as OKI were given the opportunity and the 
funding to undertake regional water quality plan- 
ning. The new 208 program (OKI 1974) differs 
from past HUD -financed water and sewer planning 
in that this new EPA -administered program deals 
with all sources of pollution, not just pollution from 
municipal sewerage systems. Other sources of 
pollution for which planning responsibility is 
given are industrial discharges and what are 
termed ir non~point ,T sources. These non-point 
sources include sedimentation and runoff of pesti- 
cides and fertilizers from agricultural areas, 
urban runoff, erosion from construction sites, 
and leachates from septic tanks. 

A common requirement of the 208 programs 
and programs of other governmental agencies 
concerned with the maintenance and control of 
water quality is the development of a knowledge of 
the interrelationships between the water quality 
parameters {turbidity, chlorophyll concentra- 
tions, etc. ) and land-use parameters (land-use 
categories and coverage, etc.)* 

To obtain this information, OKI is develop- 
ing a deterministic model capable of predicting 
water quality in rivers and lakes resulting from 
existing, and potential land-use policies. The in- 
ventory of present land use together with popula- 
tion projections will provide OKI a basis for 
developing future land-use maps. Given a future 
land use, the water quality model will be used to 
predict the imlact of future development bn 
water quality. This analysis will aid OKI in 
identifying critical areas where alternatives 
will have to be developed to. minimize any 
deleterious impact on 'water quality* This may 
involve redirecting growth to other areas where 
the impact might not be so severe, or changing 
the character of the growth to minimize any 
harmful impact on water quality* The water 
quality plan will, in effect, contain a significant 
land-us e planning element. The 208 water quality 
planning program, in fact, may provide the most 
rational basis for land-use planning available to : 
date* 
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Accurate water and land-use parameters are 
essential in the development and application of the 
water quality model. 

While many factors iaifluen.ee water quality, 
a dominating one is the use of land adjacent to and 
surrounding the lakes and rivers, the "drainage 
areas. 11 During periods of rain or thaw* this area 
discharges sediment and nutrients directly to the 
water bodies by means of surface runoff or storm 
drainage/ Each land -use category has its own 
special characteristic (EPA-143Q, 1973) which is 
important in the calculation of the quantity and 
quality of storm- water runoff . For example, 
fertilized lawns (tended grass) and paved streets 
discharge more nutrients, especially phosphorus, 
than do rangeland (untended grass) and forested 
land. Cropland is often tilled in the spring when 
rainfall is heaviest and absorbs much of the water, 
but erosion in the form of sediments, which in- . 
eludes pesticides and fertilizer, are washed into 
nearby streams. This differs from what happens 
: in a center city area where virtually all of the 
ground is covered by pavement and buildings and 
little br none of the water is absorbed into the 
earth. Instead, the water flows rapidly into storm 
sewers, carrying with it dirt from streets and 
buildings. 

To establish sediment and nutrient flows from 
the drainage areas into the waterways, accurate 
information on drainage area land use is essential. 
Land -use information presently available to plan- 
ning agencies is not adequate for water quality 
planning purposes. In almost every case, agri- 
cultural and vacant land has been lumped into one 
category such as miscellaneous. Urban land uses 
often ar e not identified (categorized) in terms 
usable for water quality planring. Also, most 208 
program planning areas are extremely large. The 
OKI Region covers 7, 024 sq km (2, 7 12 3q mi) and. 
contains 225 drainage areas. For these reasons, 
OKI decided that the traditional techniques fc-r 
land-use inventory - field inspection and interpret 
tation of aerial photographs - are impractical in 
that they ar e too costly in terms of dollar s and 
time. In its quest to evaluate new sources and 
techniques for obtaining the needed land-use infor- 
mation, OKI established and accomplished the fol- 
lowing two goals: 

* Produced a 10 -category land-use map of 
the OKI regional area showing: rangeland/ 
fallow cropland, water, cropland, two 
categorie s of for est land, cor e city / 
industrial, inner city, urban, and subur- 

. ban. The smallest detail mapped- was 
0.44 hectare (1.1 acres) and the map scale 
was 1 in. = 5, 000 ft. 
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• Produced computer tabulation of area 
covered by 16 land-use categories within 
225 dr aim ge areas and nine counties. 

OKI achieved its mapping goals within a 
period of 90 days at a cost of only $20, 000, a 
significant improvement over conventional tech- 
niques. 

OKI REGION 

The OKI objectives were achieved through 
the machine processing of the 14 April 1973 
LANDSAT-1 scene shown in Figure 1. The OKI 
region shown in this scene consists of the Ohio 
Counties of Hamilton, Clermont, Butler, and 
Warren; the Indiana Counties of Dearborn and 
Ohio; and the Kentucky Counties of Boone, Kenton, 
and Campbell. This 2, 712-sq mi region centered 
around Cincinnati is expected to increase its 
population to over 2, 000, 000 by the year 2000, a 
35% growth from its present figure of about 
1, 6GQj 000. The area is characterized by its 
many low hills and dense fores-** cover. The Ohio 
River and its tributaries drain the OKI region. 

The Ohio, Great Miami, and Little Miami rivers, 
can also be seen in the LANDS AT image of the 
region. All of the rivers and lakes in the region 
are highly valued for recreational and residential 
uses. Increasingly heavy public use makes it vital 
that water quality consideration remains as one of 
OKI's highest planning priorities. 

MACHINE PROCESSING OF LANDS AT DATA . 

The need for faster and more economical 
mapping of water quality and land use has led 
Bendbc into evaluating computer target "spectral . 
recognition" techniques as a basis for automatic 
target categorization and mapping. The cate- 
gorization techniques (Dye, 1974, 1975; Rogers, 
1974, 1975 ) have been under continued develop- 
ment at Bendix for the past 8 to 10 years, pri- 
marily using aircraft multispectral scanner data. 
More recently, LANDS AT /MSS and Skylab/EREP- 
S 192 data have been used* 

BENDIX.DATA CENTER 

The elements of the Bendbc Data Center used 
to process data for this study are shown in Fig- . 
ure 2 and include: a Bendix Datagrid® Digitizer 
System 100 for digitizing graphical data, a Bendix 
Multi spectral-Data Analysis System (M-DAS) for 
the analysis of. LANDS AT 1 'computer - compatible 
tapes" (GCTs), and a Cal Comp Plotter for the 
production of land-water categorized maps from 
the processed LANDSAT tapes. A Gerber Series 


40 Plotting Table is also used for this mapping 
function. M-DAS has been discussed (Johnson, 
1974) in detail previously. 

The nucleus of the M-DAS is a Digital Equip** 
ment Corporation PDP-11 /35 computer ■with 28K 
words of core memory, one 1. SM-word disc pack, 
two nine -track 800 bit -per -inch (bpi) tape trans- 
ports, and a DECwrxter unit* Other units are an 
Amper FR-2000 14-track tape recorder, a bit 
synchroniser and tape deskew drawers which 
can reproduce, up to 13 tape channels of multi- 
spectral data from high^ density tape recordings, 
a high-speed hard-wired special-purpose com- 
puter for processing multispectral data, a 9 1/2- 
in. drum recorder for recording imagery on 
film, and a color moving -window computer- 
refreshed display. M-DAS is the result of an evo- 
lutionary program initiated by Bendix in 1967 and 
is dedicated to the processing of remote sensing 
data. 

processing steps 

The data processing steps used (Figure 2) 
and the results achieved in transforming LANDSAT 
CCTs into the desired land-use maps and data are 
briefly summarised in the following paragraphs. 

Establish Map Categories 

The first step in the development of the OKI 
land-use map was to locate and designate to the 
computer a number of LANDS AT picture elements 
or "pixels" that best typified the land- water cate- 
gories of interest, the "training areas. " These 
areas of known characteristics were established 
from aerial photographs, and ground survey data, 
and were located on the LANDSAT CCTs by view- 
ing the tap ed data on the M- DAS TY monitor . The 
coordinates of the training areas were designated 
to the computer by placing a cursor over the de- 
sired area and assigning a training area desig- 
nation, category code, and color code. Several 4 
training areas, typically 20 to 50 pixels in size, 
were picked for each category, with, each pixel 
corresponding to a ground coverage of 57 x 79 m. 
The color code was used in later playback of the 
tap e s when the computer -categorized data are dis- 
played in the designated colors.. 

Develop Processing Coefficients 

The LANDSAT spectral measurements within 
the training area boundaries were edited by the 
computer from the CCT and processed to obtain a 
numerical descriptor (computer-processing co- 
efficients) to represent the spectral characteristics 
of each target category. ■ The descriptor s. (Dye, 
1974) included the mean signal and standard devi- 
ation for. each LANDSAT band and the c ovariance 

„„ , ORIGINAL PAGE IS 

QF POOR QUALM 


matrix taken about the mean. The desc ; iptor s 
were then used to generate a Bet of processing 
coefficients for each category. In automatic *. 
category processing, the coefficients are used by 
the computer tc form a linear combination of the 
LANDSAT measurements to produce a variable 
whose amplitude is associated with the probability 
of the unknown measurement being from the target 
sought. In category processing, the probability 
of a LANDSAT pixel arising from each one of the 
different target categories of interest is computed 
for each pixel and a decision, based on these 
computations, is reached. If all the probabilities 
are below a threshold level specified by the oper- 
ator, the computer will decide that the category 
viewed is unknown (uncategorized). 

The categorical processing of multivariate 
data at Bendix is carried out using a maximum 
likelihood procedure under multivariate normal 
hypothesis. The implementation of this procedure 
(Dye, 1975 ) for M alternatives involves evalua- 
tion of M density equations and acceptance of the 
alternative with the greatest density. The direct 
evaluation of each density equation requires (N + 

1) times N multiply and add operations if imple- 
mented with N original variables* The direct 
method, thus, is slow and costly* The procedure 
used. at Bendix greatly minimizes the computa- 
tional load with no sacrifice in computational 
accuracy* During the decision process, the 
observation is subjected to a principal component 
transformation derived fox* a given group, and 
another transformation in which the inverse of the 
covariance matrix is diagonalized with respect to 
the covariance matrix of the remaining groups. 
The first transformation gives the direction of 
maximum variance of the given group with respect 
to the background (the remaining groups). When 
the decision process is carried out, the computa- 
tion for each group now involves only N 4- 1 add 
.and. multiply operations. . In addition to the cost- 
effectiveness, it also permits an increase in 
absolute accuracy by permitting dimensionality 
reductions, since an eigen value of less than one 
from the second transformation only hinders the 
accuracy of the classification scheme and, hence, 
the associated components for an eigen value of 
less than one are rejected during computations. 

Evaluate Selection of Training Areas and 

Processing Coefficients 

Before producing categorized data for the 
entire OKI region, a number 'of tests were applied 
to evaluate the computer^ ability to perform the 
desired interpretation. The tests included gen- 
erating categorization-accuracy tables and view- 
ing the.prbcessed imagery on the M-DAS TY 
monitor. Selection of training areas, generation 
of accuracy tables, and evaluation of processing 


results through use of computer printouts and the 
TY monitor were iterative operations. The 10 
primary land -use categories listed below, in the 
order of their potential for discharging natural 
and human sources of nutrients, resulted from the 
IANDSAT processing. 

* Core City /Indus trial - Most dense indus- 
trial and commercial area, asphalt, con- 
crete, gravel, etc, , 0% vegetation* 

* Inner City - Second- most dense industrial 
and commercial area, 1 to 10% vegetation 
cover * 

• CJrban - Third-most dense industrial and 
commercial area, also includes high- 
density residential area, 10 to 50% vege- 
tation cover. 

• Suburban - Fourth-most dense area, 
housing, isolated shopping center, greater 
than 50% vegetation . cover. 

• Cropland (3 categories) . 

* Fallow Cropland. 

* Rangeland. 

• Forestland (Density 1 ) - Mixed hardwood 
(deciduous) forest and softwoods. 

• Forestland (Density 2) - Mixed hardwoods 
and softwoods, 

* Water (5 categories) - Five categories of 
water are probably; deep clear; shallow 
water; and three categories corresponding 
to different sedjtliiehfc C l icentrations,. 

In addition to the 10 major categories noted above, 
cropland was separated into three crop categories 
and water into five categories. Ground -truth acti- 
vities are underway to refine the identification of 
the crop and water categories. 

Evaluation of processed results show that 
‘some urban categories become confused with non- 
urban categories; d, e* , core city with water, 
urban with cropland, etc. However, the errors 
o c cur at : random spots as "speckles”. and, w hen 
viewed in cont ext oil a map . grid, ar e easily cor - 
recfced by an interpreter. This points up the addi- 
tional need to transfer . IANDSAT categorized data 
to a base map. 

Generate Categorised IAHD5AT Tapes 

When satisfied with the categorization accu- 
racy achieved on the land-water categories, the . 


processing coefficients were placed into the com- 
puter disk file and used to process that portion of 
the XANDSAT CCTs covering the OKI region. This 
step in the categorization processing resulted in 
new or categorized CCTs, where each LANDSAT 
pixel was represented by a code designating one of 
the 16 land -water categories. These 16 categories 
were merged into the 10 primary categories noted 
previously and mapped at a scale of 1 in# - 5, 000 
ft for the OKI region. Computer tabulations were 
also extracted from the categorized tapes to obtain 
a quantitative measure of land use within the drain- 
age areas. 

Produce Categorized Map Overlays 

To produce categorized data that will directly 
relate to a base map, the categorized CCTs were 
submitted to a second stage of processing, fa this 
stage, new tapes were generated that had data cor- 
rected for earth rotation and a format compatible 
with a computer -driven Cal- Comp Plotter, These 
tapes, when played; back by the computer, caused 
overlays of a specified land-water category to be 
drawn by the plotter (see bottom image in Figure 
2) on mylar at a scale specified by the operator. 
Examples of these mylar drawings over a map of 
drainage areas near Cincinnati are shown in Figure 
3. The examples show core city /industrial, water, 
and forestland categories at the original scale of 
1:120, 000. The overlays were photographically en- 
larged later to the final map scale of 1 in, = 5, 000 
ft, A diazo chrome material was exposed through the 
black and clear category transparencies by a litho- 
graphic plate burner and ammonia developed to pro- 
duce color-coded overlays. The color coding per- . 
mitted multiple overlays to be used simultaneously 
over the base map. 

A 10- category color coded land -use map was 
developed for the entire 7/024 sq km (2, 712 sq mi) 
region at the full JuANDSAT resolution (0* 44 
hectare or 1. 1 acre pixels ) and at a scale of 
1 in. = 5, 000 ft. The cost for this map was approxi- 
mately $2. 50 per square mile, . 

Four parameters were used in performing 
the geometric correction of LANDSAT data: space- 
craft heading, spacecraft earth latitude, adjusted 
scan line length, and spacecraft altitude. The first 
three parameters only were available on the lAKD 1 - 
SAT CCT, The heading, latitude, and adjusted 
scan line length were used to generate incremental ; 
coordinate translations of the XANDSAT data, scan 
line by scan line, to obtain along-track and cross- 
track corrections. The coordinate translation 
Increases as the computer moves through the tape,; 
and the rate of increase is a function primarily of 
spacecraft latitude. : Provisions Were made to, vary 
the coordinate translation with adjusted scan line 
length, but experience has shown that the scan line 


length varies little and the correction for scan line 
length has little effect on the data* An important 
parameter that was not on the tape is exact satel- 
lite altitude* This information is available in the 
Goddard processing facility when the tape is gen- 
erated, but for some reason it is not placed on the 
tape, Bendix has attempted to establish a channel 
through which this information can be made avail- 
able for tapes delivered to Bendix, but has been 
unsuccessful to date. Altitude variations can cause 
approximately ±0*5%. variation in the cross-track 
scale and must be corrected* The abroach used 
by Bendix was to generate one trial overlay (such 
as water boundaries) after all other corrections 
are made, overlay the trial overlay on a . map, 
select at least two control paints common to both 
the overlay and the map, and calculate a cross- 
track correction factor to compensate for satellite 
altitude variations, 

Area Measurement Tables 


LANDSAT coordinates and transformation matrices 
thus obtained are approximate, based on the use of 
the nominal spacecraft parameters* The approxi- 
mately-derived LANDSAT coordinates and transfor- 
mation are used* however, to identify the actual 
IANDSAT coordinates associated with the GCPs* 

To accomplish this, the coordinates of a GCP is 
input to the Bendix M-DAS* The approximate trans- 
formation computes the IANDSAT coordinates and 
displays the area on the TV monitor* Positional 
errors of the GCPs displayed to the operator are 
designated by a cursor to the computer, which uses 
the error measurement to derive an improved set 
of coefficients for the transformation matrix. This 
procedure is repeated on additional GCPs until the 
desired geometric accuracy is achieved. This 
rapid interactive procedure is essential for pro- 
ducing a transformation matrix that provides an 
accurate transformation of earth to IANDSAT co- 
ordinates* 


Computer -generated area measurement 
tables, illustrated in Figure 4, were produced from 
the categorized data tapes to determine land use 
within the drainage areas. To accomplish this 
step, a procedure was developed by which the 
drainage area boundaries in earth coordinates 
(latitude and longitude) are first digitized (Figure 
2) from watershed maps. The resulting digital 
tape is processed on M-DAS (Figure 2) to trans- 
form the earth coordinates to LAND SAT coordi- 
nates and to extract and tabulate land use from the 
categorized tape. The area measurement table 
provides the amount of land that falls within a. part- 
. icular category in terms of square kilometers, 
acres, and per centage of the total drainage area 
processed* These data provided a direct and use- 
ful input to OKI r s water quality model where each 
category is assigned a n loading function* 1 or pol- 
lution ecruivalent(E PA -1430, 1973). Eventually, 
by multiplying the extent of land use acreage in an 
area by these loading functions, the water quality 
model will calculate the extent of pollution emana- 
; ting from these non- point sources* 

Earth. to ERTS Coordinate Transformation - 
Theru were three basic steps involved in the auto- 
matic referencing of ground coordinates to LAND- 
SAT coordinates. The first rtep consisted ox auto- 
matic retrieval of the latitude and longitude of 
carefully selected ground control paints (GCPs) 
from a map through a digitizing process* . The. cri- 
teria for selecting these GCPs is that they can be 
easily and accurately identified on IANDSAT 
imagery. The second step consisted of converting 
the latitude and longitude of these GCPs to LAND- . . .. . 
SAT coordinates by using a theoretical transforma- 
tion derived from known and assumed spacecraft 
parameters including; heading, . s can rate, altitude, 
and a knowledge of earth rotation parameters. The 


CONCLUSION 

Machine processing of LAND SAT data pro- 
vides a rapid and economical means of mapping 
land use in watersheds of lakes and rivers. 

Although additional improvements can and are 
being made in processing techniques to increase 
mapping rates and accuracy and to reduce cost, 

OKI has demonstrated the technique s and utility of 
LANDSAT for mapping watershed land use on an 
operational basis* 

As there is still some confusion between some 
urban and non -urban categories as indicated by 
random mis clas sification (speckling), additional 
processing refinements are needed to improve 
separability of these categories. 

Machine-assisted interpretation of IANDSAT 
tapes was found to be very fast. The analysis phase 
required about one day per LANDSAT scene* Once 
the analysis was completed and the processing co- 
efficients were computed, the categorized tape was 
produced for a full IANDSAT GCT (2, 500 square 
nautical miles) in less than 30 minutes* Boundaries 
of water drainage areas were manually digitized 
.from maps of the OKI area at a rate of about six. 
per hour* The computer extracted and tabulated 
land use within these areas at a rate of one every 
3 minutes. The major time-consuming step in the 
production of the IANDSAT products was the gen- 
eration of the map overlays for the region, which 
required about 3 weeks. In the near future, equip- 
ment and techniques will be in us e that will permit 
this step to be accomplished within a day* 

The machine-processing techniques permitted 
the nine -county OKI region to be mapped to 1 -acre 
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detail for $20, 000, With application of conven- 
tional techniques, it is not uncommon for a county 
to spend this much or more to map similar cate- 
gories v/ithin a much smaller area* Additionally, 
conventional techniques based on manual interpre- 
tation of photography and field checks typically 
require a year or two to obtain a mapping product 
similar to that obtained by OKI within 90 days, 
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Develop Earth to LANDSAT 
Coordinate Transformation 

• Digitize Ground Control Points 

• Designate Location of 
Training Areas 

• Digitize Boundaries of Areas 
for which Area Printout 

Tables Are Required; Watersheds, 
Counties, Townships, etc. 



Produce LANDSAT Categorized Tapes 

• Define Land Water Categories 3nd Locate Corresponding 
Training Areas within LANDSAT Tapes. 

• Compute Category Characteristics. 

• Evaluate Training Area Selection. 

• Transform LANDSAT Tapes into New Set 
of Tapes where Each Pixel Is Coded 

to Correspond to Interpreted Land Water 
Categories. 



Generate Data and Map Products from 
LANDSAT Cc^egori zed Tapes 

• Produce Transparent Color 

Coded Overlay for Each Category; Typical 
Scales of 1 .24,000, 1 :62,500, and 1 :250,000. 

• Generate Color-Coded Imagery Where Color 
Is Used as a Code to Designate Categories. 

• Produce Tabular Computer Printouts 
Listing Area Covered by Land-Water 
Categories within Specified Political and 
Geographic Boundaries in Percent Coverage 
per Category, Acres, and Square Kilometers. 


Figure 2. Machine Processing of L-ANDSAT Da;a. 


j. 
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wmk 


Core City/Industrial Category o*er Watershed Map 


Forestland Category over Watershed Map 


Examples of Land Use Mapped from LANDSAT Categorized Tapes 
(Original Scale: 1 inch = 2 miles). 












Percent of 
Watershed 


Square 

Kilometer* 


Core City /Industrial 

Inner City 

Urban 

Suburban 

Cropland 

Fallow Cropland 

Rangeland 

Forestland 1 

Forestland 2 

Water 

Uncategorized 


Total 


Figure 4. Map of Watersheds in Southern Kenton County, Kentucky, and 

Tabular Printouts Produced from LANDSAT Categorized Tapes. 


Square 

Kilometer* 


Percent of 
Watershed 


Acre* 


Category 


1.12 


Core City /Industrial 
Inner City 


0.16 


0.0 


0.16 


58 .13 


0.24 


Urban 


1.12 


0.0 


0.16 


Suburban 


12.29 


0.05 


1.73 


Cropland 
Fallow Cropland 
Rangeland 


140.65 


0.57 


19.87 


148.68 


0.60 


20.98 


173.27 


0.70 


24.45 


Forestland 1 


168.80 


0.68 


Forestland 2 


23.82 


Water 


3.35 


0.01 


Uncategorized 


0.47 


708.75 


2.87 


100 


Total 


Category 

[percent of 
Watershed 

Acres 

Square 

Kilometers 

Core City/Indu*trial 

0.01 

1.12 

0.0 

Inner City 

0.01 

1.12 

0.0 

Urban 

6.58 

960.27 

3.89 

Suburban 

0.77 

112.91 

0.46 

Cropland 

1.59 

231.41 

0.94 

Fallow Cropland 

27.96 

4.082.55 

16.52 

Rangeland 

27.56 

4.024.42 

16.29 

Forestland 1 

13.01 

1,900.42 

7.69 

Forestland 2 

22.31 

3,258.66 

13.19 

Water 

0 

0 

0 

Uncategorized 

0.21 

30.19 

0.12 

Total 

100 

14,603.07 

59.10 






